Search for: All records

Abstract In order to better connect core-collapse supernova (CCSN) theory with its observational signatures, we have developed a simulation pipeline from the onset of the core collapse to beyond shock breakout from the stellar envelope. Using this framework, we present a 3D simulation study from 5 s to over 5 days following the evolution of a 17M_⊙progenitor, exploding with ∼10⁵¹erg of energy and ∼0.1M_⊙of⁵⁶Ni ejecta. The early explosion is highly asymmetric, expanding most prominently along the southern hemisphere. This early asymmetry is preserved to shock breakout, ∼1 day later. Breakout itself evinces strong angle-dependence, with as much as 1 day delay in the shock breakout by direction. The nickel ejecta closely tail the forward shock, with velocities at the breakout as high as ∼7000 km s⁻¹. A delayed reverse shock forming at the H/He interface on hour timescales leads to the formation of Rayleigh–Taylor instabilities, fast-moving nickel bullets, and almost complete mixing of the metal core into the hydrogen envelope. For the first time, we illustrate the angle-dependent emergent broadband and bolometric light curves from simulations evolved in 3D in entirety, continuing through hydrodynamic shock breakout from a CCSN model of a massive stellar progenitor evolved with detailed, late-time neutrino microphysics and transport. Our case study of a single progenitor underscores that 3D simulations generically produce the cornucopia of observed asymmetries and features in CCSNe observations, while establishing the methodology to study this problem in breadth. 

Abstract The interaction of supernova ejecta with a surrounding circumstellar medium (CSM) generates a strong shock, which can convert ejecta kinetic energy into observable radiation. Given the diversity of potential CSM structures (arising from diverse mass-loss processes such as late-stage stellar outbursts, binary interaction, and winds), the resulting transients can display a wide range of light-curve morphologies. We provide a framework for classifying the transients arising from interaction with a spherical CSM shell. The light curves are decomposed into five consecutive phases, starting from the onset of interaction and extending through shock breakout and subsequent shock cooling. The relative prominence of each phase in the light curve is determined by two dimensionless quantities representing the CSM-to-ejecta mass ratioη, and the breakout parameterξ. These two parameters define four light-curve morphology classes, where each class is characterized by the location of the shock breakout and the degree of deceleration as the shock sweeps up the CSM. We compile analytic scaling relations connecting the luminosity and duration of each light-curve phase to the physical parameters. We then run a grid of radiation hydrodynamics simulations for a wide range of ejecta and CSM parameters to numerically explore the landscape of interaction light curves, and to calibrate and confirm the analytic scalings. We connect our theoretical framework to several case studies of observed transients, highlighting the relevance in explaining slow-rising and superluminous supernovae, fast blue optical transients, and double-peaked light curves. 

Abstract The ongoing LIGO–Virgo–KAGRA observing run O4 provides an opportunity to discover new multimessenger events, including binary neutron star (BNS) mergers such as GW170817 and the highly anticipated first detection of a multimessenger black hole–neutron star (BH–NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH–NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH–NS mergers. Adopting optimal binary properties, a mass ratio ofq= 2, and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary’s evolution from premerger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of ∼−15 a few hours after the merger, longer than previous estimates, which did not consider the first principles–based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of ≳1 Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within ∼200 Mpc. The BH–NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals. 

Abstract Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they become core-collapse supernovae (SNe). The resulting dense circumstellar medium may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomenon were limited to 1D, leaving the 3D convective RSG structure unaddressed. UsingFLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward-moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2–4 km s⁻¹convective motions seed the intrinsic convective instability associated with the high-entropy material deep in the envelope, enabling gas from deep within the envelope to escape and increasing the amount of ejected mass compared to an initially “quiescent” envelope. The 3D models reveal a rich density structure, with column densities varying by ≈10× along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition. 

Abstract We present the first numerical simulations that track the evolution of a black hole–neutron star (BH–NS) merger from premerger tor≳ 10¹¹cm. The disk that forms after a merger of mass ratioq= 2 ejects massive disk winds (3–5 × 10⁻²M_⊙). We introduce various postmerger magnetic configurations and find that initial poloidal fields lead to jet launching shortly after the merger. The jet maintains a constant power due to the constancy of the large-scale BH magnetic flux until the disk becomes magnetically arrested (MAD), where the jet power falls off asL_j∼t⁻². All jets inevitably exhibit either excessive luminosity due to rapid MAD activation when the accretion rate is high or excessive duration due to delayed MAD activation compared to typical short gamma-ray bursts (sGRBs). This provides a natural explanation for long sGRBs such as GRB 211211A but also raises a fundamental challenge to our understanding of jet formation in binary mergers. One possible implication is the necessity of higher binary mass ratios or moderate BH spins to launch typical sGRB jets. For postmerger disks with a toroidal magnetic field, dynamo processes delay jet launching such that the jets break out of the disk winds after several seconds. We show for the first time that sGRB jets with initial magnetizationσ₀> 100 retain significant magnetization (σ≫ 1) atr> 10¹⁰cm, emphasizing the importance of magnetic processes in the prompt emission. The jet–wind interaction leads to a power-law angular energy distribution by inflating an energetic cocoon whose emission is studied in a companion paper. 

Abstract Identifying the sites of r-process nucleosynthesis, a primary mechanism of heavy element production, is a key goal of astrophysics. The discovery of the brightest gamma-ray burst (GRB) to date, GRB 221009A, presented an opportunity to spectroscopically test the idea that r-process elements are produced following the collapse of rapidly rotating massive stars. Here we present James Webb Space Telescope observations of GRB 221009A obtained +168 and +170 rest-frame days after the gamma-ray trigger, and demonstrate that they are well described by a SN 1998bw-like supernova (SN) and power-law afterglow, with no evidence for a component from r-process emission. The SN, with a nickel mass of approximately 0.09 M_⊙, is only slightly fainter than the brightness of SN 1998bw at this phase, which indicates that the SN is not an unusual GRB-SN. This demonstrates that the GRB and SN mechanisms are decoupled and that highly energetic GRBs are not likely to produce significant quantities of r-process material, which leaves open the question of whether explosions of massive stars are key sources of r-process elements. Moreover, the host galaxy of GRB 221009A has a very low metallicity of approximately 0.12 Z_⊙and strong H₂emission at the explosion site, which is consistent with recent star formation, hinting that environmental factors are responsible for its extreme energetics. 

Abstract Study of the double-detonation Type Ia supernova scenario, in which a helium-shell detonation triggers a carbon-core detonation in a sub-Chandrasekhar-mass white dwarf (WD), has experienced a resurgence in the past decade. New evolutionary scenarios and a better understanding of which nuclear reactions are essential have allowed for successful explosions in WDs with much thinner helium shells than in the original, decades-old incarnation of the double-detonation scenario. In this paper, we present the first suite of light curves and spectra from multidimensional radiative transfer calculations of thin-shell double-detonation models, exploring a range of WD and helium-shell masses. We find broad agreement with the observed light curves and spectra of nonpeculiar Type Ia supernovae, from subluminous to overluminous subtypes, providing evidence that double detonations of sub-Chandrasekhar-mass WDs produce the bulk of observed Type Ia supernovae. Some discrepancies in spectral velocities and colors persist, but these may be brought into agreement by future calculations that include more accurate initial conditions and radiation transport physics. 

ABSTRACT A subset of Type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies (∼1052 erg) that cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) that suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae form a rapidly rotating and strongly magnetized protoneutron star (PNS), a proto-magnetar. Jets from such magnetars can provide the high kinetic energies observed in SNe Ic-bl and also provide the connection to GRBs. In this work, we use the jetted outflow produced in a 3D general relativistic magnetohydrodynamic CCSN simulation from a consistently formed proto-magnetar as the central engine for full-star explosion simulations. We extract a range of central engine parameters and find that the extracted engine energy is in the range of 6.231 × 1051−1.725 × 1052 erg, the engine time-scale in the range of 0.479−1.159 s and the engine half-opening angle in the range of ∼9°−19°. Using these as central engines, we perform 2D special relativistic (SR) hydrodynamic (HD) and radiation transfer simulations to calculate the corresponding light curves and spectra. We find that these central engine parameters successfully produce SNe Ic-bl that demonstrates that jets from proto-magnetars can be viable engines for SNe Ic-bl. We also find that only the central engines with smaller opening angles (∼11°) form a GRB implying that GRB formation is likely associated with narrower jet outflows and Ic-bl’s without GRBs may be associated with wider outflows. 

ABSTRACT Neutron star merger accretion discs can launch neutron-rich winds of >10−2M⊙. This ejecta is a prime site for r-process nucleosynthesis, which will produce a range of radioactive heavy nuclei. The decay of these nuclei releases enough energy to accelerate portions of the wind by ∼0.1c. Here, we investigate the effect of r-process heating on the dynamical evolution of disc winds. We extract the wind from a 3D general relativistic magnetohydrodynamic simulation of a disc from a post-merger system. This is used to create inner boundary conditions for 2D hydrodynamic simulations that continue the original 3D simulation. We perform two such simulations: one that includes the r-process heating, and another one that does not. We follow the hydrodynamic simulations until the winds reach homology (60 s). Using time-dependent multifrequency multidimensional Monte Carlo radiation transport simulations, we then calculate the kilonova light curves from the winds with and without dynamical r-process heating. We find that the r-process heating can substantially alter the velocity distribution of the wind, shifting the mass-weighted median velocity from 0.06c to 0.12c. The inclusion of the dynamical r-process heating makes the light curve brighter and bluer at $$\sim 1\, \mathrm{d}$$ post-merger. However, the high-velocity tail of the ejecta distribution and the early ($$\lesssim 1\, \mathrm{d}$$) light curves are largely unaffected.